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Abstract 

Quadratic operators are used in transforming the model Hamiltonian H of one correlated and 
dispersive band in an unique positive semidefinite form coopting both the kinetic and interacting 
part of the H . The expression is used in deducing exact ground states which are minimum energy 
eigenstates only of the full Hamiltonian. It is shown in this frame that at half filling, also dispersive 
bands can provide ferromagnetism in exact terms by correlation effects . 



I. INTRODUCTION 



Exact results for non-integrable many body models open new perspectives in the de- 
scription and understanding of physical systems of large interest, or physical situations 
capturing attention by unusual properties. Metallic ferromagnetism 1 , superconductivity in 
ladder systems^, exotic phases in polymer chains^, subtle conducting and insulating phases 
in rare earth compounds 4 , spontaneous ferroelectricity^, stripes or checkerboards in two 
dimensions 6 , or metal-insulator transition in two dimensional disordered systems^ are some 
examples in this respect. 

On the technical side this field is in fact almost completely open, hence several procedures 
have been tested or are under development. Without to be exhaustive, we exemplify below 
some of these, as for example process definition followed by unapproximated treatment of 
the model conditions^, Hilbert space properties captured in the low density limit 9 , equal 
rigorous upper and lower bounds for the ground state energy^, unitary transformation of 
eigenfunctions and Hamiltonians 1 ^, Schrieffer- Wolff transformations effectuated up to infi- 
nite order—, the optimal ground state method 1 ^, or the use of positive semidefinite operators 
and related techniques^ 4 -. In the present paper we concentrate on, and develop the last of 
these methods since the procedure based on positive semidefinite operators often penetrates 
in other techniques as well. 

By using positive semidefinite operators (P) in deducing exact ground states (\^ g )), 
one first casts the Hamiltonian (H) in a positive semidefinite form H = P + C, where 
C is a constant depending on the coupling constants of the starting H. Since zero is the 
minimum possible eigenvalue of P, the unapproximated ground state is obtained from the 
requirement P\^ g ) = 0. In the process of transforming the starting Hamiltonian in a 
positive semidefinite form, mostly one finds or uses P = Y^T ax Pn, where the number n max 
of different types of positive semidefinite terms often can reach even n max = 8 — 10, their 
structure being diverse, containing linear, bi-linear, cubic or quartic combinations of the 
starting fermionic operators^. The requirement for the ground state in the presence of 
different P n contributions becomes P n \^ g ) = for all n G [l,n max ], hence the deduced \^ g ) 
represents separately the ground state of different parts of the Hamiltonian, often being the 
ground state separately of the kinetic and interacting part of H, which reduces considerably 
the application possibilities of the method. 
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In this paper we overcome this inconvenience by using quadratic, e.g. nonlinear combi- 
nations of the starting fermionic operators in defining new non-fermionic operators A\. Here 
the index i represents a lattice site, but Ai is not local, has a finite extension since captures 
contributions also from the neighborhood of the site i. The operators A± are used in the 
construction of an unique G = J2i A\A^ positive semidefinite form, leading to the expression 
of the transformed Hamiltonian H = G + C, containing an unique positive semidefinite oper- 
ator, namely G (the contribution C, as before, being a constant). Excepting the test of this 
procedure taken at infinitely large on-site Coulomb repulsion^, for finite and nonzero value 
of the interaction this approach is for the first time applied here in this paper. The gain one 
achieves is that now an unique positive semidefinite operator concentrates both the kinetic 
and interaction part of H. Hence the corresponding ground states are not independently 
eigenstates of the kinetic and interacting part, or independently eigenstates of different parts 
and terms from H, but only eigenstates with minimal energy of the full Hamiltonian. 

Relating the studied Hamiltonian, one concentrates here on one-band correlated itinerant 
systems holding as well spin-spin interactions. The model Hamiltonian is in this case of 
extended Hubbard type (denoted hereafter by Hq) containing also spin-spin interactions 
Hs, hence H = H + Hg holds. A such Hamiltonian structure is particularly often debated 
in connection to Kondo lattices 17 or t — J type of models- , but generally emerges in the 
study of strong correlation effects when additional interactions between nearest neighbours 
are taken into consideration^, or real materials^Saii are described. We further note that 
in the presence of strong correlations, similar models have been also used in the study of 
charge-ordered states, and metal-insulator transitions in non-ferromagnetic cases^. 

The here deduced exact eigenstates characterize dispersive bands at half filling and are 
localized and fully saturated ferromagnetic ground states only of the full H. At this point we 
remark that the exact ferromagnetic ground states at finite concentration of carriers deduced 
up to this moment are or related to a specific band structure (flat band ferromagnetism^, or 
Lieb's degenerate middle band ferromagnetism in a bipartite lattice^), or appear given by 
the ability of the interaction to select the ground state from the existing minimum energy 
eigenstates of the kinetic (not necessarily one-particle^) parlr^. Our results show that 
arbitrary dispersive half filled bands can be ferromagnetic in exact terms. Since only the 
whole Hamiltonian provides these ground states, these not emerge given by a specific band 
structure, nor appear as interaction selected ground states from existent eigenstates of other 
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parts of the Hamiltonian, but are created by both kinetic and interaction contributions, 
hence their eigenstate nature disappears when the interactions are turned off, or when the 
kinetic part is neglected (e.g. in the localized limit). 

The remaining part of the paper is structured as follows. Section II. describes the model 
used, Section III. presents the exact transformation of the Hamiltonian in positive semidefi- 
nite form, Section IV. gives the general expression of the deduced ground states and presents 
the proof of the ground state nature, Section V. describes particular cases when the deduced 
ground states hold and analyzes the physical properties of the solutions. Finally, Section 
VI. concluding the paper, closes the presentation. 

II. THE USED HAMILTONIAN 

Our starting Hamiltonian H = H + Hg describes in one dimension a correlated band 
characterized by 

Ho = ^(t a cl a c i+1 ^ + H.c.) + U^n^n- hl + Y^ ^1,02^,0^+1,^, 

i,cr i i o-i,(72 

+ ^(^l,a^i-a4,aCi+l,a + ^2,a^i+l- CT c|+l, CT Ci, CT + ^-C.), (1) 

i,cr 

where c\ a creates an electron with spin a at the site i, n^ a = c\ a c- lta represents the particle 
number operator, t a is the hopping matrix element, U represents the Hubbard term, V aii(T2 
gives the density- density nearest-neighbor interaction strengths, while Wi t<7 , 1 = 1,2 describes 
the correlated hopping. 

Inside the band, also a nearest neighbor spin-spin interaction is present given by 




where 5 ,± = S x ±iS y , Si = J2 a ,/3 A,oPo-,^,Pi an< ^ &a,p describes the Pauli matrices. During 
this paper periodic boundary conditions are used. 

III. THE TRANSFORMATION OF THE HAMILTONIAN 

In order to deduce exact ground states for the Hamiltonian of the problem one transforms 
H into an unique positive semidefinite form. In order to do this we introduce a bi-linear 
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operator A ; defined by 

+ WiCj+l.tCi,! + W 2 Ci+l,lCi,-\ + W3Ci+l, T Ci + i,| + W4Ci,TCi,|, (3) 

where the pref actors z 9 ,wg, 9 = 1, 2, 3, 4 are considered site independent numerical parame- 
ters. One observes that J2i A\A{ reproduces the operators contained in Eqs.( fTf2l . and indeed 
one finds that the exact mapping 

# = ^iJii-aiiV-ao (4) 

i 

holds, if the numerical coefficients zg, wg, a a , are defined via the following matching con- 
ditions 

4 

A 0=1 

t+ j_ + + 

Vt,t + x = n,4 + ^ = 0, 

t / I 1 2 1 I 1 2 t / ^ I 1 2 1 1 1 2 

v u - "j = N + Kl , vi,t - -7- = N + |Wi| , 

1 i2ii 1 2 1 1 i2ii i2 ^-L ^ ^ 

(7 = \Z 3 \ +\Z 4 \ +\W 3 \ + |W 4 | , — = W X W 2 + ^1 ) 

+ 2*2:3 = W 2 ,i, W3W2 + Z2Z3 = —W 2 ,'[, 

Wlw 2 + -22-24 = —Wl f l, wlWi + Z{Z 4 = W\^, 

W3W4 + z\z 3 = F, w{z 2 + w* 2 zx + wlz A + w\z 3 = 0, (5) 

where represents the number of lattice sites, and N is the total particle number operator. 
Since the total number of particles is a constant of motion, via Eq. (jlj) we matched the starting 
Hamiltonian into an unique positive semidefinite form using the conditions presented in 
Eq.©. 

IV. THE EXACT GROUND STATE AT HALF FILLING 

Let us consider the bond operator 

B\ = «t4t + a l d ll + PA+1,1 + PA+i,l> ( 6 ) 
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where a a and /3 a are numerical prefactors. One observes that {B}, Bj} = holds for all i, j 
and all a a , (3 a . 

Let us analyze the wave function 

N A 

i^) = nA t io), (7) 

i=l 

where |0) is the bare vacuum with no fermions present. Since in Eq.([7j) one introduces N\ 
electrons into the system, \ty g ) is defined at half filling. 
If for all i one has 

M*g) = 0, (8) 

then clearly, \^ g ) represents the ground state wave function of H. Since the fermionic 
operators from A\ in Eq.Q are present in Y\- } Bj only in Bj, Bj ±1 , it results that in order to 
satisfy Eq.Q one must has for all i sites 

AbUbLbIM = 0, (9) 

from where one finds 

ZxPi - z 3 ai = 0, + z 3 ai; = 0, z x a^ - z 4 f3^ = 0, z 2 a L + z±P± = 0, 

wi(api - a^a^) + w 3 (/3fa± - PiP^) = 0, 

w 2 (a 2 fr - a^Pi) - w 3 ((3fa^ - /^ftaj) = 0, 

w 4 (a^ T - a^aiP{) + Wx(f3 2 a^ - A/ 5 !"!) = °> 

-w 4 (a 2 p l - a T a^ T ) + w 2 (f3 2 a i - /3 T /3^o; T ) = 0, 

-w 4 (/3jar T - PiPiati) - WiPlfr - w 2 0\P\ = 0, 

w 4 ((3 2 a l - /3 T /3|Q; T ) - w 1 [3 2 [3 l - w 2 pp l = 0, 

wz(a 2 \P^ — a^a^Pi) — wia 2 a-\ — w 2 a 2 a^ = 0, 

— w 3 (a 2 Pi — a-^aiPi) — Wia 2 a^ — w 2 a 2 cti = 0. (10) 

In Eq. lflO]) the first row comes from the coefficients of the terms containing three creation 
operators created by the product A\B\, while the remaining equations are provided by the 
rest of the terms from Eq.([9]). 

Based on Eqs. pilBIh the ground state energy is given by E g = — CLq — CLi N. 
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V. THE DEDUCED SOLUTIONS 



In the following one presents the obtained solutions for the system of equations Eqs. (l5][T0]) 
denoted by Solutions I. and II. which describe two different regions of the phase diagram. 

A. Solution I. 

Taking into consideration W\, a = W 2 , a — W, the parameters of the Ai, B\ operators 
become 

ia Z l -ia Z l ( \q\ 2 ~ 1 \ 

z 2 = -z x , z 3 = qzie , z A = — e , wi = — ( 



q 2 q 

^lQ —in in 

w 2 = -iwi, w 3 = e , w 4 = z x e , 

q 

«T = Xa h A = ^a|e JQ , (3 l = qa^ 10 , (11) 

where z±, a, a b A, q remain arbitrary coefficients at the level of ( TTTT) . The Hamiltonian pa- 
rameters turn to have the property t = t a , V\ = V a>a , V 2 = V a - a , a = UN^, ai = —U, and 
one has 

Xir 1 T/ J * T/ U+J * T U no\ 

W = ~2' V i = ~> V2 = ^~' J± = ~2~' ( ^ 

where U > 0, t and J z can be arbitrary chosen, and the parameters |zi|,|g|,a becomes 
determined by \zi\ = \t\/y/U, (1 + |<?| 2 )/|<?| = |t|/|zi| 2 , and e~ ia = (q/\q\)sign(t). The 
coefficients A and «j remain arbitrary, and the operator constructing the ground state wave 
function from ([7]) becomes 

B\ = ai [(cl + Acf )T ) + qe ia (ct +u + Acf +ljT )]. (13) 

B. Solution II. 

If one considers W\ t(T = Wi, W 2ta = W 2 and takes v = z 3 /w A = —w^/zl non-zero and 
real, the parameters of Ai, B- t operators become 



z 2 = -zi, z 3 = Pe ia zi, z A = —p- z i> Wi = -zi, 

-ia i 

w 2 = zi, w 3 = -v——zi, w A = -Pe zi, 

P V 

ai = \a b fr = \Pe* a ai , f3 i = Pe ia a b (14) 



where P = y\v\/ (y/2 + sign(v)), a is given by sign{t) = —e~ ia , and z±, a;, A, v are arbitrary 
non-zero coefficients. Furthermore one has a\ = a /N\ = ^W\W2/2(2 + P 2 + P~ 2 ), and the 
Hamiltonian parameters become 

Wi = -sign{t y^ VP, W 2 = -sz#n(t)V2>i| 2 \u\VP^, 
J ± = -pWyW*, V 1 = V 2 = ^ - J ± , 

V2 V k| 1 1 ; |z/|V2WW 2 M 

where i/, J 2 remain arbitrary, sign(Wi) = sign{W2) must hold, and one has \zi\ 2 = 
^W\W2/2. The operator present in the expression of the ground state wave function (j7j) is 
given by 

B\ = a^Xcl + cfj - Psign(t)(\c\ +1A + cJ +1 J], (16) 
so its form is similar to that given in Eq. ({T3l) . 



C. Physical properties of the deduced solutions 

For the ground state (J7J), in both obtained cases (I13|16p describing different regions of 
the phase diagram, the B\ operator is constructed from l\ = (Acj ^ + cj ^) type of blocks, for 
which Jj Jj = holds for arbitrary A. Consequently, apart from a normalization constant, 
n£i4 + |o) from © becomes Ili=a -m|0)j hence the normalized ground state in both cases of 
Solutions I. and II. can be written as 

N A 
i=l 

where A is arbitrary, and iV = in (fTTj) fixes the filling at half. Denoting by (...) = 
(ty g \...\ty g ) ground state expectation values, and using as usual S? = 1/2(^1 — nyj, S3~ = 
(S~)\ = Cj+Ci^, = ± iSf, and S" = X^Si for the spin operators, one finds (S 2 ) = 
(N\/2)(N\/2 + 1). This result is motivated by the fact that since A is the same on each i, 
the local spin orientations are the same on each site. Hence the ground state represents a fully 
saturated ferromagnetic state. Furthermore, since only local Hamiltonian terms contribute 
in the ground state energy, and independent on i, (n^) = 1 — (^4) = |A| 2 /(1 + |A| 2 ), 
(^ijT^i,;) = 0; the ground state is localized 4 . Consequently, the obtained ground state at 
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half filling of the studied one (dispersive) band model is a localized and fully polarized 
ferro magnet. 

Note the important aspect that \^f g ) is the ground state only of the full Hamiltonian. 
Hence not a specific band structure produces here the ferromagnetic behaviour (as for ex- 
ample in the case of the flat band ferromagnetism^, or Lieb's degenerate middle band 
magnetism in a bipartite lattice^), nor the interactions select the ground state from the 
minimum energy eigenstates of the kinetic part^^. In the described case the ferromag- 
netism at half filling can appear in exact terms in an arbitrary dispersive band being only 
the ground state of the full Hamiltonian, e.g. is provided by an incontestable correlation 
effect. 

As described previously in details^, for a non-integrable model as that studied here, exact 
ground states can be obtained only on different cuts of the phase diagram corresponding 
to different decompositions in positive semidefinite form of the starting Hamiltonian. Since 
for the presented two solutions in Sections V.A, V.B, the expression of the operator Ai is 
different, these solutions correspond to different decompositions, hence provide the deduced 
ground state in different regions of the phase diagram. The fact that the same solution 
appears in different regions underlines the stability of the ferromagnetic phase, since finite 
H parameter modifications lead back to the same ordered phase. However, the ferromagnetic 
localized phase has a finite extension. This is seen from the fact that since \^f g ) is not an 
eigenstate of the non-interacting Hamiltonian, turning off the interactions one recovers the 
conducting behaviour of a non-interacting dispersive band, e.g. a paramagnetic metal to 
ferromagnetic insulator transition must be present in the phase diagram of the system. 

We further note that resembling ferromagnetic states have been found in similar systems 
by other methods as wel&-^^-. 

VI. SUMMARY AND CONCLUSIONS 

Single band extended Hubbard models containing spin-spin interactions are studied by 
a positive semidefinite operator technique, leading to ferromagnetic and localized exact 
ground states in different regions of the phase diagram via quadratic operators used in the 
positive semidefinite decomposition leading to an unique operator term in H. Hence the 
deduced ground states are minimum energy eigenstate only of the full Hamiltonian, e.g. 
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not are connected to a specific band structure, nor are separately eigenstates of the kinetic 
and interacting parts. The ferromagnetism emerges at half filling in an arbitrary dispersive 
band, being provided by clear correlation effects. 
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